clear all; clc; close all;

%Внутренние параметры
R=220; 
C=33e-12;
L=100e-6;




omega0=1/sqrt(L*C);
f0=omega0/2/pi;
fprintf('f0=%fMHz\n',f0/1e6);


%Ось времени
T=1/(1000*f0);
t=0:T:(15*1/f0);
lt=length(t);

%Выделение памяти
Uc=nan(1,lt);
%Ul=nan(1,lt);
dUc=nan(1,lt);
i=nan(1,lt);
di=nan(1,lt);

%Внешние параметры 

E = 10*ones(1, lt);


%Начальное состояние 
Uc(1)=0;
i(1)=0;
di(1)=0;
Ul(1)=di(1)*L;



for k = 2:lt
%Ul(k) = E(k);
Ul(k) = Uc(k-1) + i(k-1)*R;
% Напряжение на катушке
di(k) = Ul(k) / L;
% вызывает приращение тока,
i(k) = i(k-1) + di(k-1)*T;
dUc(k) = i(k) / C;
% который определяет заряд конденсатора
Uc(k) = Uc(k-1) + dUc(k)*T;
end 


U = Uc;


figure(1);
plot(t*1e9, [Uc; E; U]);
xlabel('t, ns');
ylabel('U_c, E, U, Volt');
legend('U(t)', 'E(t)', 'U(t)');
grid on;





















